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ABSTRACT 

We present diffraction limited, 10/im imaging polarimetry data for the cen- 
tral regions of the archetypal Seyfert AGN, NGC1068. The position angle of 
polarization is consistent with three dominant polarizing mechanisms. We iden- 
tify three distinct regions of polarization: (a) north of the nucleus, arising from 
aligned dust in the NLR, (b) south, east and west of the nucleus, consistent with 
dust being channeled toward the central engine and (c) a central minimum of 
polarization consistent with a compact (<22pc) torus. These observations pro- 
vide continuity between the geometrically and optically thick torus and the host 
galaxy's nuclear environments. These images represent the first published mid- 
IR polarimetry from an 8-m class telescope and illustrate the potential of such 
observations. 
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Introduction 



The unified tlieory of Seyfert (Sy) type active galactic nuclei (AGN) liolds tliat all 
types of Sy AGN are essentially the same object, viewed from different lines of sight (LOS). 
Surrounding the central engine is a geometrically and optically thick, dusty, molecular torus, 
obscuring the broad emission line region from some LOS. In this scheme, the Sy classification 
depends solely on the LOS and exact torus properties. Such theories received a major boost 
through the d etection of scattered, and h ence polarized, broad emission lines in the spectrum 
of NGC1068 (lAntonucci fc Millen Il985l ). revealing an obscured Sy 1 central engine in the 
previously classified Sy 2 AGN, entirely consistent with unified theories. 

Whilst fundamental to unified theories, the torus remains difficult to image directly 
at optical/IR wavelengths, with perhaps the most direct ob servation of the torus made 
by speckle interferometry in the near-IR (jWeigelt et al.l 120041 ). Strong evidence for signifi- 
cant amounts of obscuring material in the central lOOpc-scale nuclear regions of NGC1068, 
possibly in the forna of a torus, is prov i ded by observations o f CO and HCN emission 
( jPlanesas et al.l Il99ll : iJackson et al.lll993l : ISchinnerer et al.l 12000) , and recent Chandra X- 
ray observations (e.g. lOgle et al.l l2003l ). Mid-IR sp ectroscopy reveals a moderately deep 
'l^jj^ 0.4) silicate absorption feature at the nucleus (IRoche et al.lll984j : ISiebenmorgen et al. 



20041 ), whose strength is approximately const ant up to ~1" south of the br i ghtest mid-IR 
point jMason_eia,l.ll2006l: iRhee fc Larkirul2006l ). Applying the lNenkova et al.l (120021 ) clumpy 



torus model, 



Mason et al. J2006h suggeste d the torus is corn pact (< 15pc), in good agreement 
with mid-IR interferometric observations (IJaffe et al.ll2004f). Further evidenc e for a compact 



torus was found through AO fed H2 l-OS(l) observations flDavies et al.ll2006l ). finding a 15pc 
clump of H2 extending from the nucleus at the same PA as the line of masers. The observed 
extent of the torus, or nuclear obscuring ma terial in NGC1068, is partly dependent on the 
wavelength and/or observational technique. I Young et al.l (119961 ) used imaging polarimetry 
to observe the silhouette of obscuring material against the southern ionization cone, which 
they attributed to the torus, with a derived diameter of ~200pc in the H-band. 

The close proximity (1" = 72pc) and high brightness of NGC1068 makes it an ideal target 
for polarimetry , a tra di tionally photon-starve d ap plication. Near-IR stu dies of NGC1068 by 
Packham et al.l (119971 ). iLumsden et al.l (119991 ) and lSimpson et al.l (120021 ) clearly detected the 
bi-conical ionization structure in scattered light. In the nuclear regions, there is a trend to 
a position angle (PA) of polarization being perpendicular to the radio jet with increasing 
wavelength. Modeling of the nuclear regions requires both an extended area of scattering 
particles as well as dichroic absorpti o n of nuclear emis s ion, p ossibly by dust in, or associated 
with, the torus (I Young et al.l (119951 ). IWatanabe et al.l (120031 )). 



Bailey et al.l (119881 ) found that the PA of polarization rotates by ~70° between 4 and 
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5/im. The lO/iin spectropolarimetry of lAitken et al.l (119841 ) showed a similar PA of polariza- 
tion to that at 5/im, and a constant degree of polarization through the silicate absorption 
feature. These results are entirely consistent with the predicted 90° change from dichroic 
absorption to dichroic emission from aligned dust grains. That the PA change was only ~70° 
is attributable to dilution of the dichroic emission component by polarized flux in the ex- 
tended scattering c ones, most likely fr om dichroic emission from dust in the narrow emission 
line region (NLR) JSailev et al.lll988h . 



To investigate the contri butions of the v arious polarizing mechanisms and structures 
in the nucleus of NGC1068, iLumsden et al.l (119991 ) performed imaging polarimetry using 
a broad-band 8-13 fim filter. These data represented the first and only published mid-IR 
imaging polarimetry of an AGN, but their interpretation was complicated by the ~0.7" 
resolution of the data. To take advantage of the improved spatial resolution attainable from 
an 8 m-class telescope, we obtained new mid-IR imaging polarimetry of NGC1068 during 
commissioning of this mode at the Gemini North telescope. 



2. Observations 

We obtained imaging polarimetr y of NGC1068 during commissioning of the polarimetry 



unit of Michelle ( IGlasse et al. 1119971 ) on UT 2005 December 19 on the Gemini North 8.1m 
telescope. These observations were primarily aimed at measuring the degree and position 
angle (PA) of polarization with NGC 1068 as a test object, and hence used a limited on-source 
time of 148 seconds. Michelle uses a Raytheon 320 x 240 pixel Si:As IBC array, providing 
a plate scale of 0.1" per pixel in imaging mode. Images were obtained in the 9.7/im {SX 
= 1.0/im, 50% cut-on/off) filter only, using the standard chop-nod technique to remove 
time- variable sky background, telescope thermal emission and so-called "l/f detector noise. 
The chop throw was 15" and the telescope was nodded every ~90 s. The chop throw was 
fixed at 0° (N-S). Conditions were photometric and the observations were diffraction limited 
(~0.30"FWHM). 

Michelle employs a warm, rotatable half wave retarder (or half wave plate, HWP) to 
modulate the polarization signal, located upstream of the entrance window of the dewar. A 
cold wire grid polarizer is used as the polarimetric analyzer, located in a collimated beam. 
Images were taken at four HWP PAs in the following sequence: 0°, 45°, 45°, 0°, 22.5°, 67.5°, 
67.5°, 22.5° in the first nod position, and the sequence repeated in the second nod position. 
In this manner, the Stokes parameters can be computed as close in time as possible, reducing 
the effects of variations in sky transmission and emission. This sequence, however, requires 
many motions of the HWP, and is therefore under evaluation with a view to reducing the 
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number of HWP motions to increase observing efficiency. Data were r educed using the 



Gem ini IRAF package in conjunction with Starhnk POLPACK software (IBerry fc Gledhil] 



20031 ). The difference for each chop pair in a given nod position and HWP PA was calculated, 
and then differenced with the second nod position at the same HWP PA. Images were aligned 
through shifting by fractional pixel values to account for slight image drift between frames, 
and then the Stokes parameters I, Q and U computed using POLPACK. A total of 20 nod 
positions were recorded, and residual array /electronic noise was removed through use of a 
median-filter noise mask. The data were reduced through creation of four individual I, Q 
and U maps and also through coadding all frames at their respective HWP PA first and 
then producing a single I, Q and U map. The S/N in the latter method is slightly higher in 
the individual Q and U maps, presumably due to a 'smoothing' of the Q and U during the 
co-addition; these are the data used in this paper. 

The efficiency and zero angle calibration of the polarimeter were measured through 
observations of tw o polarized sources and comparison with measurements published by 



Smith et al.l (120001 ). The instrumental polarization was estimated to be <0.3% through 
observations of two stars that fulfilled the criteria of (a) high proper motion (hence nearby), 
(b) high galactic latitude (to minimize the presence of Galactic dust) and (c) an intermediate 
spectral type star (to minimize intrinsic stellar nebulosity). 



3. Results 

Figure 1 shows the total fiux image (color-scale and contours) with the polarization 
vectors overlaid. The polarization vectors are plotted where the S/N is >54 in total fiux, 
and contours are linearly spaced in intensity, starting at a S/N of 27. Figure 2 shows the 
polarized intensity map, produced by multiplying the degree of polarization by the total 
intensity image. As in Figure 1, only where the S/N in the total intensity image is >54 
are data plotted. The resultant polarization vectors are contained within an approximate 
N-S oriented ellipse, major/minor axes 1.7"/1.2" respectively. The integrated degree of 
polarization within that ellipse is 2. 48 ±0.5 7% at a PA of 26.7±15.3°. The errors in the 
degree and PA of polarization are estimated through independent measurements of the four 
individual polarization maps and computing the standard deviation. It should be noted that 
the exposure time in those four individual maps was very short, where systematic effects could 
dominate, and hence the quoted error s may be an ove r estim ation. The degree of polarization 



is higher than the 1.30% measured by lLumsden et al.l (119991 ) in a 2" aperture, consistent with 



an increased observed polarization as often arises with improved spatial resolution. The PA of 
polarization is significantly different from Lumsden's measurement of 49° in a 2" aperture. 
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However, our data shows a PA rotation of 94° from the K„ band data of Packham et al. 



(119971 ) and lLumsden et aLl (119991 ). entirely as expected if the dominant polarizing mechanism 
changes from dichroic absorption to emission between the two wavebands, as described in 
311 We speculate the iLumsden et al.l (119991 ) ~0.68" results suffered significantly greater 
contamination in their beam, possibly from surrounding exten ded polarization , as co mpared 
to our ~0.30" results. Additionally, the wider bandwidth of iLumsden et al.l (119991 ) would 
have been more affected by the different and competing polarizing components. 

The degree and PA of polarization suggests contributions from three components. The 
first extends ~1" north of the mid-IR peak and is coincident with the inner regions of the 
radio jet, with a PA of polarization approximately N-S. The second region extends south, 
east and west of the nucleus, with a PA of polarization of ~35°. Finally, the degree of 
polarization drops to a minimum very close to the mid-IR total flux peak, which we believe 
is most likely to arise from an unresolved polarization contribution with a PA of polarization 
approximately orthogonal to that of the more extended emission, leading to a reduction in the 
measured polarization. The polarized intensity image reveals polarized emission extending 
north of the mid-IR peak, and two areas of enhanced flux east and west of the mid-IR peak, 
and a minimum close to the mid-IR total flux peak. Table 1 summarizes the locations and 
polarization components. 



Discussion 



Polarization at mid-IR wavelengths most likely arises from either dichroic absorption or 
emission, both due to dust grains with a preferred alignment. The integrated PA of polar- 
ization in these data confirms and enhances the interpretation of the near-90° PA flip from 
near- to mid-IR wavelengths, with t he factor ^2.5 i r icrease in sp atial resolution providing a 
more accurate and consistent result. iGalliano et al.l (l2003l . 120051 ) suggested, based on spatial 
coincidence, the [OIII] clouds in the ionization cone are the dominant mid-IR sources away 
from the compact torus. The polarized flux image shows a similar spatial correspondence 
with the [QUI] , and the PA of po larization north of the nucleus is consistent with the in- 
terpretation of iBailey et al.l (119881 ) of dichroic emission in the NLR, possibly through dust 
alignment via jet streaming or a helical magnetic field associated with the jet. We discount 
directly observed synchrotron radiation from the radio jet accounting for the polarization, 
as an extrapolation of the radio emission to the mid-IR provides too little flux. Hence, this 
data confirms the extended mid-IR polarized emission north of the nucleus is dominantly 
from dust in the ionization cone. 



South, east and west of the nucleus, as the PA of polarization is perpendicular to that 
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in the near-IR where the polarization is thought to be produced by dichroic absorption, we 
suggest the dominant polarizing mechanism is dichroic emission by grains aligne d to the same 



field d irection as the absorbing g rains, in agreern ent with other authors (i.e. iBailey et al. 



( 119881 )). Due to the Barnett effect (ILazarianl (120031 ). and references therein), grains align with 
their short axes parallel to the local magnetic field, and the PA of polarization is parallel 
to the direction of the magnetic field for dichroic absorption and orthogonal for dichroic 
emission. The location of the polarized emission areas and PA of polarization is suggestive 
of warm aligned dust grains being channeled from the host galaxy towa rd the torus. I ndeed , 
the PA of the polarized fiux is coincident with the H2 material that iDavies et al.l (120061 ) 
associated with molecular material in a compact torus. 

Dichroic absorption in an unresolved optically thick central region could account for 
the minimum in polarization close to the peak of mid-IR flux, with a PA of polarization 
approximately orthogonal to the more extended dichroic emission to the east, west and 
south. Alternatively the mid-IR flux in the innermost regions could arise from a strong 
mid-IR, intrinsically unpolarized, source. However, there is tentative evidence of the central 
regions showing a twist in the PA of polarizatio n, tending toward a sim ilar PA found in the 
dichroic absorption at near-IR wavelengths (e.g. iPackham et al.l (119971 )). which supports the 
dichroic absorption interpretation. In both possibilities, a compact (<0.3" (<22pc) diameter) 
torus could account for this result. If correct, the polarization minimum indicates the true 
position of the central engine, which is not coincident with the mid-IR total flux peak, but 
displaced by ~0.2" to the west. 



CO (ISchinnerer et al.l I2OOOI ) and optical HST (ICatchpole fc Boksenberd 119971 ) obser- 
vations are interpreted as evidence of a warped molecular disk on lOOpc scales, partially 
obscuring the nuclear regions o f the host galaxy and ionization cone pointing away from 
Earth. Indeed, ISchinnerer et al.l (120001) speculate this material, ra ther than a compact torus, 
is responsible for obscuring the AGN. Elitzur fc Shlosman ( 20061) suggest the lOOpc molec - 
ular structure is the e x tension of the pc sca le disk of masers ( IGreenhill &: Gwinru ( 119971 ) , 
Gallimore et al.l (l200ll ). iGalliano et al.l (120031 )). We suggest that our data provide continuity 
between the geometrically thick torus (height /radius ~1) to the flatter, larger galactic disk 
(height /radius ~0.15). 

The western polarized featur e is considerably l arger than the compact (< few pc) torus 



Jaffe et al. 


2004; 


Mason et al. 


2006 



f 

Ip 



Radomski et al.ll2006l ) on the bas is of mid-IR i r naging and modeling, but much smaller than 
the suggested torus detected by lYoung et al.l (119961 ). However, the feature is detected in 
polarized flux, a technique that increases contrast by removing the dominant, unpolarized, 
emission. Distinct from total flux, polarimetric observations are therefore potentially much 
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more sensitive to emission from the putative faint, outer regions of the torus where the 
interaction with the inner regions of the host galaxy must occur. We suggest that a way 
to reconcile the evidence for a compact torus w i th these observatio ns and others, such as 



extended sihcate absorption (IRoche et al.l (120061 ) . iRoche et al.l (120071 )) and 100 pc-scale CO 



discs, is that the compact, geometrically and optically thick torus is often surrounded by a 
larger, more diffuse structure, associated with the dusty central regions of the host galaxy. 
Where the torus ends and the host galaxy dust structure starts may be more of a question of 
semantics rather than a true physical boundary. These observations examine the interaction 
between the host galaxy and possible entrainment into the outer torus regions. Further 
multiple-wavelength polarimetric observations of both NGC1068 and other AGN are required 
to test this hypothesis. 

We are grateful to the Gemini, UKIRT and ATC science and engineering staff for their 
outstanding work on Michelle and the Gemini telescope, and wish to note especially Chris 
Carter. Based on observations obtained at the Gemini Observatory, which is operated by 
the Association of Universities for Research in Astronomy, Inc., under a cooperative agree- 
ment with the NSF on behalf of the Gemini partnership: the National Science Foundation 
(United States), the Particle Physics and Astronomy Research Council (United Kingdom), 
the National Research Council (Canada), CONICYT (Chile), the Australian Research Coun- 
cil (Austraha), CNPq (Brazil) and CONICET (Argentina). 
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Table 1: Polarization component summary 



Locale, Aperture, A 


Degree of Polarization 


PA of Polarization 


Emitting Component 


Nucleus, 1.7"xl.2" 


2.48±0.57% 


26.7±15.3° 


Several 


9.7/im 








North region 


-2% 




Dichroic emission from dust 


9.7/im 






aligned through jet interaction 


East, West, South 


-3.5% 


-35° 


Dichroic emission from 


regions, 9.7/im 






galactic dust or torus outer edge 


Innermost region 


<1% 




Dichroic absorption or 


9.7/im 






unpolarized source 


Nucleus, 2" 


1.3±0.05% 


49±3° 


Several, including dichroic 


10/im 






emission 


Nucleus, 2" 


4.11±0.46% 


120.6±2.38° 


Several, including dichroic 


2.2/im 






absorption 



°9.7/xm data from this paper 
^10/im data from lLumsden et al. 
^2.2/im data from iPackham et al. 



(Il999l) 
(|l997t ) 



Fig. 1. — Total flux image (color) with the polarization vectors for the central regions of 
NGC1068. The length of the vector is proportional to the degree of polarization, and the 
angle shows the PA of polarization. Each pixel is 0.1", and the 10% polarization scale bar is 
shown in the upper right. North is up, and east is to the left. 



- 12 - 




